Dissipative quantum dynamics of bosonic atoms in a shallow ID optical lattice 
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We theoretically study the dipolar motion of bosonic atoms in a very shallow, strongly confined 
ID optical lattice using the parameters of the recent experiment [Fertig et al, Phys. Rev. Lett. 
94, 220402 (2005)]. We find that, due to momentum uncertainty, a small, but non-negligible, atom 
population occupies the unstable velocity region of the corresponding classical dynamics, resulting 
in the observed dissipative atom transport. This population is generated even in a static vapor, due 
to quantum fluctuations which are enhanced by the lattice and the confinement, and is not notably 
affected by the motion of atoms or finite temperature. 
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The rich phenomenology of the dynamics of an atomic 
Bose-Einstein condensate (BEC) in a periodic optical lat- 
tice potential has recently attracted considerable theo- 
retical and experimental interest. In particular, the non- 
linear mean-field interaction of the BEC may give rise 
to dynamical and energetic instabilities in the transport 

properties of the atoms HEUSEHHHHE^- M ° St 
experiments to date have been carried out in shallow lat- 
tices where the motion of the BEC can be accurately 
described by the classical mean- field theories. In the 
strongly interacting regime, with very high lattice poten- 
tials, the atoms were shown to undergo a quantum phase 
transition to the Mott insulator state [ill IT^| . Recently 
strongly damped oscillations of bosonic atoms were ex- 
perimentally observed by Fertig et al. at NIST in a 
tight elongated ID atom trap, even in the limit of a very 
shallow lattice potential. In this Letter we numerically 
investigate the damped oscillations of bosonic atoms in 
a shallow lattice using the parameters of Ref. and 
show that the observed damping rate can be explained 
by the growing population in dynamically unstable ve- 
locity region, resulting in phase decoherence. This pop- 
ulation appears even without the dipolar motion and is 
not notably affected in the studied case by the atom dy- 
namics. Moreover, for fixed nonlinearity, both, the pop- 
ulation at unstable velocity and the damping rate, have 
approximately exponential dependence on the atom num- 
ber. For larger atom numbers the quantum fluctuations 
are suppressed and the finite-temperature effects can be- 
come observable. 

In the recent NIST experiment 01 an array of inde- 
pendent ID atom 'tubes' were generated by applying 
a strong transverse 2D optical lattice potential. The 
bosonic atoms in each tube were also trapped along the 
axial direction by a harmonic magnetic potential and by 
a very shallow periodic optical lattice. The dipole oscil- 
lations of atoms along the weak axial lattice were excited 
by suddenly displacing the harmonic trap by means of 



applying a linear magnetic field gradient. In the absence 
of the lattice, the dipolar motion is undamped and un- 
affected by the interactions. However, significantly in- 
hibited dynamics was observed even in the presence of 
a very shallow lattice with the height of 0.25E r , where 
E r = 2n 2 h 2 /m\ 2 is the lattice photon recoil energy with 
the wavelength A. The strong effect of such a shallow 
lattice on dissipation is surprising, since the lattice mod- 
ulates the atom density only by 6%. Due to the limited 
resolution of the imaging system, it was not possible to 
measure the damping of the oscillations in the center-of- 
mass (cm) position of the atoms, but the oscillations in 
the cm velocity were imaged after a time-of- flight, follow- 
ing the sudden turn-off of the trapping potentials. 

Here we numerically study the dipole oscillations of 
bosonic atoms in a ID lattice within the truncated 
Wigner approximation (TWA) for the complete multi- 
mode field operator, beyond the discrete tight-binding 
approximation. The previous theoretical studies of the 
quantum dynamics of bosonic atoms in a combined har- 
monic and lattice potential indicate that quantum fluc- 
tuations can result in strong dissipation [bj, |ljj |ljj, |lj| • 
In Ref. ^3 the dynamics was studied within the tight- 
binding approximation to the TWA for deep lattices. 
This approach is much less accurate for very shallow lat- 
tices studied here and in Ref. Moreover, we find that 
the tight-binding approximation to the TWA for shallow 
lattices strongly underestimates the generated noise and 
the effect of vacuum fluctuations on the dynamics. 

In Refs. [l3L H§| the atoms were confined in a 2D ar- 
ray of decoupled tight ID tubes by means adiabatically 
loading the atoms into a strong 2D optical lattice. We 
write the atom density in a ID tube at (xq, j/o) : 
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where p 2 = (x—xo) 2 + {y—yo) 2 , Ri are the Thomas-Fermi 
radii of the 3D BEC, with R = {R x R y R z ) 1 / 3 , N a denotes 
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the total atom number, and d — A/2 the lattice spac- 
ing. By expanding around the lattice site minimum we 
may estimate l± ~ (S/mui) 1 '' 2 by the harmonic trap fre- 
quency in the transverse direction lo±_ = 2y / ~sE r /h, where 
s denotes the lattice height in the units of E r . In Ref . , 
s = 30 and A = 810nm, yielding uj±_ ~ 2tt x 38kHz. For 
N a = 1.4 x 10 5 atoms 13], R x ~ 14/xm, i? a ~ 
i2 2 ~ 10.6/im, and within the classical radius there are 
TTR x R y /d 2 ~ 5400 atom tubes. From Eq. (JJJ we obtain 
the number of atoms in each tube N(xo,yo) = J d 3 rp = 
5N d 2 (l - xfjRl - y 2 /R 2 y f/ 2 /2iiR x R y . At the central 
tube, iV(0, 0) ~ 65 atoms. 

When a shallow optical lattice is applied, the wave 
function Eq. (JIJ is modified along the z axis. In the 
experiment the displacement of the harmonic potential 
was about S ~ 3/im~ 2.21, where I = (hjrniS) 1 ! 2 and the 
trap frequency u> ~ 2tt x 60Hz. Within 2R Z , there are 
about 52 lattice sites. 

We study the NIST experiment of dipolar oscillations 
within the TWA, by suddenly displacing the harmonic 
trap. In the TWA, in a tight elongated trap uj <C uj±, 
the dynamics of the ID classical stochastic field iftw{ z , t) 
follows from the nonlinear equation |19|: 

i^w=(-^ + V + g\^ w \^ w , (2) 

which coincides with the Gross-Pitaevskii equation 
(GPE). Here the potential is a combined harmonic 
trap and a periodic optical lattice V(x) = muj 2 z 2 /2 + 
sE r sin 2 (nz/d) and g — 2huj±a, where a denotes the scat- 
tering length. The thermal and quantum fluctuations are 
included in the initial state of ipw i n Eq- 0) which repre- 
sents an ensemble of Wigner distributed wave functions. 
Our basic formalism of the TWA is similar to Ref. . 
Initially, before displacing the trap, the gas is assumed to 
be in thermal equilibrium and we approximate the field 
operator ip(z,t = 0), within the Bogoliubov theory [20| : 

4>(z) = V>o(z)ao + [ u j( z ) & 3 ~ v j( z ) & \] ■ ( 3 ) 

j>0 

Here ipo is the ground state solution of the GPE, ob- 
tained by evolving the GPE in imaginary time, and 
(djdo) = Nq, the number of ground state atoms. The 
mode functions Uj(z) and Vj(z) (j > 0), and the ener- 
gies 6j for excited states are obtained by solving the cor- 
responding Bogoliubov equations (ljj in the harmonic 
trap, modulated by the shallow lattice potential. In 
the TWA the time evolution of the ensemble of Wigner 
distributed wave functions [Eq. (J2J] is unraveled into 
stochastic trajectories, where the initial state of each 
realization for tpw is generated according to Eq. © 
with the quantum operators (a,-, at) replaced by complex 
Gaussian-distributed random variables (a:j , a* ) . These 
are obtained by sampling the corresponding Wigner dis- 
tributions [l5| . In particular, for the quasiparticles the 



Wigner function is that of the ideal harmonic oscillators 
in a thermal bath 21], with the width fij + \, where 

fij = (aj&j) = [exp {tj/ksT) — 1] _1 and the 1/2 width 
at T = for each mode represents the quantum noise. 

Because the Wigner function returns symmetrically 
ordered expectation values, (a*atj)w — % + | and 
(ctj)w = (oij)w = (a 2 )w = 0, for j > 0, it is gener- 
ally not possible to extract from the TWA simulations 
the correlation functions for the full multi-mode field op- 
erator [lflj . Consequently, we derive the desired normally 
ordered expectation values by defining the grou nd state 
operators cij for each individual lattice site j [l5| : a,j (t) 
j dztpQ(z)tpw(z 7 t), where the integration is over the jth 
lattice site, ipw(z,t) is the stochastic field, determined 
by Eq. (J2J), and tpo(z) is the ground state wave function, 
obtained by evolving the GPE in imaginary time. This 
provides us with the basis for ip(z) = • aj<f)j(z), where 
the ground state functions 4>j are restricted in the jth 
site, and the normally ordered expectation values are eas- 
ily obtained with respect to ip, e.g., for the cm position 
z cm , the position fluctuations Az = \{z 2 ) — (z) 2 ] 1 / 2 , and 
for the normalized phase coherence between the central 
well and its ith neighbor d = | (ajfli) \/ ^/ntflH- 

We numerically study the dipole oscillations using the 
experimental parameters |l3l |. In the central atom tube 
N ~ 65, resulting in the nonlinearity Ng ~ i2Qhwl for 
a ~ 5.313nm. In Fig. \I\we use the same fixed Ng, but 
vary N, and show the cm position z cm of the atoms for 
s = 0.25. The atoms are initially in thermal equilibrium 
at t = and at uito = 1 we displace the center of the 
harmonic trap from z = Otoz = <5 = 2.16Z~ 7 Ad. In 
the initial state Q, we generate the noise for 300 modes, 
corresponding to 4-5 lowest energy bands. Even though 
the atoms mainly remain in the lowest band, synthesizing 
the TWA noise for only 70 modes (~one band) can un- 
derestimate the damping and the phase decoherence sev- 
eral tens of percents. This emphasizes the importance 
of the multi-mode approach to the TWA, beyond the 
tight-binding approximation, and is consistent with our 
previous observations [l^. In ID TWA, unlike in 3D |22| . 
we did not find any divergence of the calculated results, 
when the mode number was increased. For each run, we 
typically use 600 stochastic realizations. The integration 
of the nonlinear evolution in Eq. @ is obtained using the 
split-step method [2^] on a spatial grid of 4096 points. 
The numerics becomes much faster for larger atom num- 
bers, due to smaller quantum fluctuations. 

As in Ref. ^j^i we model the cm motion as a damped 
harmonic oscillator z cm = —kz/m* — 2^z, so that for un- 
derdamped motion z cm (t') = — e -7 * (5 cos ftt' + c sin fit'), 
with c = 7<5/fi and il = ^/k/m* — 7 2 , where t 1 = t — to 
and m* denotes the effective mass. At t' = we have 
Zcm = —5 with v C m = 0. Even in a very shallow lattice, 
s = 0.25, we find a significant damping for N — 75 and 90 
atoms, with j/w ~ 0.020 and 0.013 [Fig.[T]. The damp- 
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FIG. 1: The oscillations of the cm position of atoms from the 
harmonic trap center at T = with fixed Ng (top left; curves 
from top N = 400,125,75), the damping rates (top right) 
7(iV) and 7- u (A r ) (upper line), and the normalized quasimo- 
mentum distributions at ujt = 8 (bottom left) with shaded 
unstable region; the inset shows the population n(p)h/l at 
v — 9uil where the classical motion becomes unstable. The 
trap is instantaneously displaced at bit = 1. The phase coher- 
ence between the initially central well and its fifth neighbor 
C 5 (bottom right; curves from top N = 1000, 400, 125, 90, 75). 

ing is clearly correlated to the loss of phase coherence. 
The dependence of 7 on N is approximately exponen- 
tial j(N) ~ O.llwexp {-N/43), resulting in 7 ~ 0.025w 
at the central tube with N ~ 65, which is close to the 
experimentally observed value of 7 ~ 0.03-0.04(±0.01)u 
[13j . We also checked that the damping rate j v , obtained 
from v cm — dz cm /dt, is close to 7. For s — 0.25 the single- 
particle effective mass is almost equal to m and we find 
fi ~ u>, with very little effects of fluctuations in fi. 

In the experiment the dynamics was measured by 
imaging the expanding atom cloud, indicating that 7 did 
not represent the value in the central tube, but the aver- 
age value over all the tubes. Although in the experiment 
g is constant in each tube (as N varies), we can obtain 
an upper bound limit for the average damping rate by 
extrapolating our fitted j(N) (for fixed Ng) to small N, 
so that for j/oj « 1, 7 = l{N)N{x l ,x j )/N a . Here 
N(xi,Xj) is the atom number in a tube at (xi,yj) and 
the summation is over all the tubes. By replacing the 
sum by an integral, we obtain 7 ~ 0.047w. 

We also studied thermal effects on the dipolar motion 
[Fig. |2| . In the experiments 0] it was not possible to de- 
rive information about temperature and the effect of ther- 
mal atoms on 7. For small N we found thermal fluctua- 
tions to be negligible compared to quantum fluctuations. 
Only for ground state atom numbers clearly exceeding 
400, we can recognize any notable finite-T effects on 7. 
This strongly suggests that the experimentally observed 
damping |l3J duly was due to quantum fluctuations. 

The classical nonlinear evolution of a BEC in a peri- 
odic lattice is susceptible to dynamical instabilities aris- 
ing from both constant velocity, larger than a critical 
value HSHI3, anc ^ acce l era ti° n or force, smaller than 




FIG. 2: Thermal effects on the coherence and the cm motion. 
The coherence C5, as in Fig.0(left); curves from top represent 
(N ,k B T/huj) = (1000,0), (900,21.3), (400,0), and (360,10.5). 
The cm (right) for the same cases. Only the finite-T, No = 
900 case has notable damping (7 ~ 0.002lj). 

a critical value @, H • The onset of the instability is 
associated with inhibited transport, fragmentation of the 
density profile, and the significant broadening of quasi- 
momentum distribution. In Ref. |l3l| the importance of 
the dynamical instabilities on dissipation was not clear, 
since the velocity was much smaller than the critical value 
and no momentum broadening was observed, despite the 
strongly damped motion. 

In the combined harmonic and lattice potential, the 
separation of the acceleration and the velocity insta- 
bilities is not straightforward, since the trap induces a 
local force muj 2 z and, e.g., applying a constant force 
moj 2 F on atoms is the same as displacing the trap: 
z 2 /2 — zF = (z — F) 2 /2 + const. Moreover, a weakly 
accelerating lattice does not equal to applying a constant 
force, since the acceleration in a shallow lattice (s = 0.25) 
does not induce notable cm motion. We numerically in- 
tegrated the classical GPE for a BEC (without quantum 
or thermal noise) and found the onset of the dynamical 
instability at 5 C ~ 8-9Z, corresponding to a maximum ve- 
locity v£ m — 8-9ul. For smaller displacements, no insta- 
bilities were observed on the time scale of the experiment 
H3, e.g., due to the weak force induced by the trap. With 
quantum fluctuations, the classically observed sharp on- 
set of the dynamical instability is smeared out due to the 
phase uncertainty of atoms between neighboring lattice 
sites 0|- This is similar to a BEC in a double- well po- 
tential, where the macroscopic self-trapping of the atom 
population can decay due to dissipation |24j . 

In order to investigate the role of the classical dynami- 
cal instabilities in the dissipative quantum dynamics, we 
calculated the quasimomentum distribution in the TWA 
from ip and found atoms occupying the dynamically un- 
stable velocity region of the corresponding classical sys- 
tem; Fig. H This population is larger and grows more 
rapidly, the larger the observed damping: e.g., at v = 9luI 
and Lut = 8, n(p)H/l ~ 0.02 and 0.08, for N = 225 
and 75, and both n(p) and 7 show similar exponential 
dependence on N. For (N,Ng/hujl) = (90,450) (here 
7 ~ 0.026w), the ratio of the population at the unstable 
velocity v ~ 9tol to that at the peak of the distribution 
was 0.07 at uit = 3 and as high as 0.09 at cot = 6. It is 
this finite high velocity occupation that results in dissi- 
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FIG 
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3: The coherence Cs (left) in a displaced trap, in a lat- 



tice moving with a constant velocity 2.16/cj, in an accelerated 
(0.27/oj 2 and 0.54Zoj 2 ) lattice, and in a static lattice. Only the 
0.54/a; 2 accelerated lattice exhibits a lower C5. The quasimo- 
mentum distribution (right) in a displaced trap (off-centered 
curve) and in a static potential at cut — 6. 



pation. Alternatively, we may say that the position un- 
certainty Az due to vacuum fluctuations is comparable 
to the critical displacement 8 C , resulting in dissipation, 
even though 8 < 8 C . In the example case Az ~ 5-6Z . 
This is consistent with the observation in Ref. that 
the damping arises from the large depletion, due to the 
tight transverse confinement, which leads to the popula- 
tion of high-momentum states. The tails of the momen- 
tum distribution can become very lon g w hen the atoms 
reach the ID Tonks- Girardeau regime |25| . 

Moreover, we also studied the atom response to an 
accelerated lattice, to a lattice moving with a constant 
velocity, and to a static lattice, without displacing the 
trap; see Fig. We used the parameters of Fig. ^ with 
N = 125 and Ng = 405. The constant velocity was 
v = 2.16lcu and the acceleration 0.27lcu 2 over the time 
period of tot = 8, so that the final velocity was 2.16lcu. 
In all the cases z CTn remained close to zero, without no- 
table cm motion. The decrease in d in all the three 
systems was almost the same as for the displaced trap, 
indicating that the motion of neither the atoms, nor the 
lattice, had significant effect on the phase coherence in 
the experiments. The occupation of the dynamically un- 
stable region in the quasimomentum distribution in the 
static potential was approximately equal to that in the 
displaced trap case (e.g., at p = ±9h/l, n ~ 0.052/fi), re- 
sulting in a similar number of high velocity atoms in the 
two cases. Moreover, the variation of this atom number 
was very little affected by the dipole oscillations. Con- 
sequently, neither phase decoherence nor the population 
with unstable velocities, which result in the damping of 
the dipole oscillations, are here due to the motion of the 
atoms, but the combined effect of the strong transverse 
confinement and the lattice. This is contrary to the sys- 
tem studied in Ref. 

Hi 

where the atom motion itself was 
argued to generate the decoherence. 

As expected, both the increase of the lattice height and 
the initial displacement enhance dissipation. We studied 
the dynamics for s — 1, 8 = 2.16?, and obtained for 
(N,Ng/fkjl) = (330, 350), (430, 340), (830, 330), 7/w ~ 
0.019,0.012,0.003, at T = 0. Also with s = 1 the damp- 
ing is notable for 8 much smaller than the classical critical 



displacement S c ~ 6-71. On the other hand, for N = 175 
and Ng ~ 380So;/, the increase of 8/1 from 2.16 to 4.32 
increased j/u from 0.004 to 0.009. 

In summary, we studied dissipative atom transport in 
a shallow ID lattice. At a fixed nonlinearity, smaller TV 
results in larger phase and momentum uncertainty and, 
consequently, larger population in the dynamically un- 
stable velocity region. Surprisingly, this population is 
generated due to the confinement and the lattice, even 
with neither the atoms nor the lattice moving. The qual- 
itative agreement between the TWA and the experiment 
[lflj . despite the approximations involved in generating 
the initial state at small N, also represents a consider- 
able success of the TWA which is traditionally consid- 
ered more suitable in the limit of large field amplitude 
[rjj.l22T| . For instance, the time-dependent Hartree-Fock- 
Bogoliubov theory predicts much too small a damping 
[16j . This suggests the TWA could provide a powerful 
technique to study vacuum fluctuations also in several 
other related condensed matter and atomic systems, such 
as corrugated superconducting nanowires 26]. 
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After our work was completed, the dynamics in a deep 
lattice was analyzed in the fermionized limit |27j . in the 
Mott state. In this limit the TWA with a bosonic field 
amplitude is no longer valid. 
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